Quasideterministic realization of a universal quantum gate in a single scattering process 
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We show that a flying particle, such as an electron or a photon, scattering along a one-dimensional waveguide 
from a pair of static spin- 1/2 centers, such as quantum dots, can implement a CZ gate (universal for quantum 
computation) between them. This occurs quasi-deterministically in a single scattering event, hence with no need 
for any post-selection or iteration, and without demanding the flying particle to bear any internal spin. We show 
that an easily matched hard-wall boundary condition along with the elastic nature of the process are key to such 
performances. 

PACS numbers: 03.67.Lx, 03.67.Hk, 42.50.-p 

Interfacing static qubits ^ mediated by flying particles is 
a prominent paradigm in the quest for efficient ways to im- 
plement quantum information processing (QIP) [2, 3]. As a 
major motivation, this is the only way to jointly address quan- 
tum registers located far from each other, thus featuring no 
direct mutual interaction (this is usually sought to favor lo- 
cal addressing). Within this general framework, over the past 
few years a research line has thrived around the idea that the 
crosstalk between the static qubits can be mediated by parti- 
cles scattering from them ll4l- [l5ll . 

Yet, all of such strategies unavoidably face an inherent ma- 
jor drawback. For a given quantum task to be efficiently ac- 
complished, the link between the static objects should occur 
by means of the local interaction of each static object with a 
quantum flying bus. Namely, this should feature inherently 
quantum motional and/or internal (pseudo) spin degrees of 
freedom (DsOF). To do so, however, the coupling between 
the flying and static particles will in general entangle them so 
as to bring about decoherence affecting the DsOF of the static 
objects. Owing to such effect, the attainment of satisfactory 
figures of merit thus demands for further actions to comple- 
ment the above interaction. These typically comprise iterated 
injection of the flying particles and post-measurements over 
their DsOF lul- USll . While such bus-dynamics conditioning 
generally enhances the performances, it usually comes at the 
cost of making the process probabilistic and may be demand- 
ing in practice (e.g. spin post-selection of mobile-electrons in 
semiconducting media lllql ). 

Moreover, as far as scattering scenarios are concerned, 
there appear intrinsic hindrances to accomplish certain tasks 
such as the implementation of a two-qubit quantum gate 
(TQG), typically the most challenging and essential process in 
most QIP architectures especially when allowing for universal 
quantum computation (QC). To see this, assume that we need 
to realize a TQG between a flying qubit and a scattering center 
(SC ) endowed with spin in a one-dimensional (ID) waveguide 
II17I1 . Even if the dynamics is conditioned to either the reflec- 
tion or transmission channel, the resulting process lacks uni- 
tarity, i.e. a paramount prerequisite for a TQG, unless quite 
specific regimes of parameters and, importantly, interaction 
models are addressed lUTIl . Analogous considerations a for- 



tiori hold when many scattering centers are present and a gate 
involving their DsOF is sought, which will be our focus in this 
work. Scattering-based scenarios thus appear as adverse are- 
nas to perform quantum algorithms, which arguably explains 
why mere entanglem ent g eneration was almost exclusively in- 
vestigated to date iRl- USll . Nevertheless, scattering -based im- 
plementations are attractive because of the low demand for 
control. One normally just needs to set the itinerant-bus wave 
vector and wait for the collision to occur, thus bypassing any 
interaction-time tuning (usually a significant noise source). 
Further benefits such as the resilience to relevant detrimen- 
tal factors including static disorder II13I1. ph ase noise Ill4ll and 
imperfect particle-wave-vector setting II15I1 have been shown. 
Except for the attempt in Ref. iflTIl such advantages have so 
far been harnessed solely for mere entanglement generation 

Here, we discover a simple strategy for the realization of 
quantum gates between static qubits through a particle scat- 
tering from them. The injection of the latter, which is not 
demanded to bear any internal DOE, followed by its multi- 
ple scattering from the SCs suffice to qnasi-dete rministically 
achieve the gate in one shot. Also, neither post-selection of 
any kind nor repeated sending of the flying mediators are 
required. Thus, besides unveiling an unexpected suitability 
of scattering-based methods i4 r-15-1 to achieve unitary oper- 
ations, our work sets a milestone within the distributed-QIP 
context l2,|3tl- For two static qubits a universal CZ gate natu- 
rally arises showing the effectiveness of our scheme. 

Central idea. Consider a monochromatic spinless particle 
/ of wave vector k propagating along a ID wire that impinges 
on an array of SCs [see Fig. 1(a)]. Once multiple scattering 
has occurred and assuming, importantly, that the process is 
elastic, /can only be found either reflected or transmitted with 
wave vectors-^ and k, respectively [see Fig. 2(b)]. Let {|v)) be 
a basis of the SCs' Hilbert space and \^) one of its elements, 
whereas \+k) are momentum eigenstates of/. Let !*!',„) = \k)\ii) 
be the overall system's initial state. As / is scattered off the 
final state reads 
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FIG. 1. (Color online) Scheme working principle, (a) / impinges on 
a set of SCs in p. (b) After the interaction, / is either reflected or 
transmitted while the SCs undergo a non-unitary quantum map. (c) 
With a perfect mirror beyond the centers, / can be only back reflected 
and the unitary R is applied to the SCs. 



where ry^ (fy/j) is a reflection (transmission) probability am- 
plitude corresponding to the initial and final centers' states 
l/i) and |v), respectively. Defining a reflection (transmission) 
operator R (T) in the Hilbert space of the SCs such that 
<v|^|;u) = Ry^ = ry^i «yu|f |v) = Ty^ = fy^) Eq. ([B can be ar- 
ranged as l^'f) = \-k}R\jj} + \k}T\iJ.}. Tracing over/, the final 
SCs density operator reads .R|yL(){ju|.R' +7'|/i)(/i|7''. Thus when 
the centers are initially in an arbitrary (in general mixed) state 
p (in general mixed) their final state is given by 



p' = RpR"'' + TpT 



(2) 



The normalization condition 2y(ky^P + |fy;jP) = 1 Vp entails 
RR'' + TT^ = 1, where 1 is the identity operator of the SCs' 
Hilbert space. We would like the scattering process to imple- 
ment a multi-qubit gate, which is unitary, between the SCs. In 
general, this is not the case as is evident from Eq. ^ showing 
that the SCs undergo instead a quantum map |1] comprising 
reflection and transmission channels. 

Such hindrance can be overcome in a very natural fashion 
by simply inserting beyond the centers a perfect mirror [see 
Fig. 1(c)]. As this introduces a hard- wall boundary condition 
(BC) preventing / from trespassing the right end, this way the 
transmission channel is in fact fully suppressed. Hence, fy^ = 
and (|2]i reduces to 

p = RpR"'" , 

where now RR' — WR - 1, i.e. in the presence of the mir- 
ror R becomes a unitary gate. Thereby, for a monochromatic 
wavepacket such gate is deterministically implemented when- 
ever / scatters from the centers. Remarkably, this holds re- 
gardless of the specific scattering potential. Rather, this affects 
only the type of achieved gate. 

Having illustrated how a hard-wall BC guarantees the pro- 
cess unitarity, it is now natural to wonder whether there exist 
elastic one-channel scattering processes allowing for multi- 
qubit gates universal for QC tlj]. We will show that this is 
indeed the case. With this aim, we first focus on a simple 
paradigmatic setting, setup A, comprising two spin- 1/2 cen- 
ters, i.e. two qubits lT|], each coupled to a massive particle 
embodying /. We identify a regime such that a CZ gate, uni- 
versal for QC Jl]], naturally arises. We next address setup 
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FIG. 2. (Color online) Setups implementing our scheme, (a) An 
electron along a quantum wire with DQDs. (b) A photon along 
a waveguide with A-type atom-like systems. The right end of the 
wire/waveguide behaves as a perfect mirror. 



B, comprising multi-level atom-hke systems and photons and 
within experimental reach in several scenarios Ill8l422ll . We 
show that this can work as an effective emulator of setup A. 

Setup A. This setting [see Fig. 2(a)] compiises two identical 
spin- 1/2 scatteiing centers, 1 and 2, lying along a semi-infinite 
wire along the x-axis at x = xi and x = X2, respectively, each 
coupled to a scattering particle / of mass m. The wire ends at 
X = X3. Let {|0),, |1),) be an orthonormal basis for the /th cen- 
ter (i - 1,2). In practice, one can consider a semiconducting 
quantum wire or a carbon nanotube ||23I1 where an electron 
populating the lowest subband can undergo scattering from 
two double quantum dots (DQDs) 11241] to which it is electro- 
statically coupled. As shown in Fig. 2(a), each single-electron 
DQD is in |0) (1 1)) if the upper (lower) dot is occupied, hence 
implementing an effective qubit 112411 (tunneling between the 
upper and lower dots is negligible). Also, the electrostatic 
coupling is negligible for state |1). The Hamiltonian is thus 
modeled as (we set ft = 1 throughout) 



H^^+ r6(x-xi) |0)i<0| + r6(x-X2) |0>2< 
2m 



(3) 



where p = -id/dx is the momentum operator of/, while F is 
the height of each contact potential scattering centered at x = x,- 
(/ - 1,2) |25| . Note that the scattering potential in (O is dis- 
persive because it cannot induce either 1 or 2 to flip between 
|0) and |1). Upon scattering, each initial SCs' state |q'iQ'2)i2 
(or,- = 0, 1) simply but crucially picks up its own phase shift. It 
is trivially checked that for a given state a = {ai, 02} H takes 
the effective form Ha - §;;, + 'Zii=i,2^SaiO S{x-Xi). The prob- 
lem thus reduces to a particle scattering from spin-less poten- 
tials. We label with ra the /'s reflection probability amplitude 
corresponding to Ha, where the subscript here specifies both 
the initial and final centers' state (these coincide owing to the 
dispersive interaction). Likewise, in the DQDs computational 
basis 1100), |01), |10), |ll))Rnecessarily has the diagonal form 
R = diag(roo,roi,rio,rii). 

Adopting a standard procedure, to derive ra we assume 
that / is left-incoming and seek the stationary state ^aix) 
fulfifling Ha'i'a(x) = k^ / (2m)'i' a{x) of the form '¥a{x) = 
^a+{x)+^a-{x) with (for simplicity we drop the dependance 
on a whenever unnecessary) 

'i'4x)^-^m-x)+ai[e{x)-e{x-X2)]+a2e{x-X2)}e"'\{4) 
V2^ 

'¥-{x)^-={r0{-x)+biie{x)-0{x-X2)]+b20(x-X2)}e-'''-\i5) 
V2^ 



where A: >0, 9{x) is the Heaviside step function and we set xi - 
0. Due to Eqs. (|4| and (|5]l, ^+{x) [^^(x)] represents the right- 
propagating (left-propagating) part of ^(x). Note that ^{x) is 
specified by the five or-dependent coefficients {r,a\,b\,a2,b2]. 
These are found by requiring that ^(x) and its derivative with 
respect to x ^''{x) match the five BCs 



^>{xT)^^>(x^) (/=i,2), >i'(x3) 

^'V'\,. = 2m V6„.Q ^(jc,) (/ = 1 , 2) . 



0, 



(6) 
(7) 



Eqs. (|6]l ensure the matching of the wave function at the cen- 
ters' locations (first two) and the hard-wall BC owing to the 
end of the wire at x = X3 (latter equation). Eqs. Q are stan- 
dard constraints on the discontinuity of ^'(x) at the centers' 
positions AT|.v, =*I"(x+)-*i''(x7) due to the ^-potentials in ^ 
[they are derived by integrating the Schrodinger equation (SE) 
corresponding to Ha over an infinitesimal range across x = xi 
and x = X2]. Solving the linear system (l6]l-(|7]i we end up with 

ror = - exp{2/ arg [exp(-/A:x3i)-2y5„,o[cos kx2\ 

-(i+2y6a^^ sinA:x2i]sinA:x32 -2y5„,osin^X3i]} , (8) 

where x,j = x,- - Xj and y - mT/k. As expected, this yields 
\ra\ = 1 regardless of all the parameters and a, namely / is 
reflected back with certainty. Consider the regime 



A:x2i — nn, kx22 — (n' + l/2)7r y s> 1 



(9) 



where n and n' are arbitrary integers. Replacing ( fT6] l in (|8]l 
and using that A:x3i =A:x32+A:x2i, we end up with 



'"00 = roi = rio 



-ru 



-1 



(10) 



which yields the gate matrix R = diag(l, 1,1,-1) (up to an ir- 
relevant global phase factor), i.e. the well-known CZ gate j^. 
This proves that a single scattering process can implement a 
universal TQG. Intuitively, unlike Ref. il7il here the geome- 
try secures unitarity leaving the physical parameters free to be 
tuned so as to match a CZ. 

Setup B. In this setting [see Fig. 2(b)] / is a photon prop- 
agating along a ID waveguide, having geometry analogous 
to setup A, and scattering from two three-level atom-like sys- 
tems. Each "atom" / = 1,2 has a A-type energy-level con- 
figuration consisting of a twofold-degenerate ground doublet 
spanned by states {|go), l^i)) and an excited state \e) [see 
Fig. 2(b)]. Hence, a two-photon Raman transition between 
|§o) and \g\) can occur through absorption and re-emission 
of a scattering photon. Assuming a linear photon dispersion 
law (jL>-vk with co the photon energy and v the group veloc- 
ity, the free-field Hamiltonian in the waveguide reads IuTIuSII 
Hf - -iYjch± ^ dxvdc\{x)dxCd{x) with v+ - -V- - v and 
c|(x) [ci(x)] the bosonic operator creating a right (left) prop- 
agating photon at x. The free atomic Hamiltonian reads 
Ha - Wo 2,=i_2 \^)i{^V where ojq is the energy gap between \e) 
and the ground doublet. The field-/th atom coupling is mod- 
eled as m Hfi = J ^ dxS(x-Xi)[c(x)S]+n.c.] (under the 
usual rotating-wave approximation), where c(x) = c+(x)-i-c_(x) 
annihilates a photon at x regardless of its propagation direc- 
tion, J is the rate associated with each transition |^/),- <-> |e), 
(Vy = 0, 1) and S'- - Ziy |e),(gy|. The full Hamiltonian thus 



reads H - H j+Ha+H jx+H f2- It is convenient to use as basis of 
each ground doublet the symmetric and antisymmetric com- 
binations of the two ground states \4>^)i-{\go)i±\gi)i)l^. As 
10"); is a dark state 112811 . i.e. S .\4>^)i - 0, the atomic raising 
operator takes the effective form S . = yl2\e)i{(f>^\. Thus the 
Raman process does not couple |0^) and |0"). It should be 
clear now that by taking |0) = \(p^) and \\)-\(p^) for each atom, 
as long as these are initially in the ground doublet, setup B in 
fact possesses all the key features of A. Indeed, if the /th atom 
is in 1 1), = \4>^)i the corresponding potential 77/, vanishes. If in 
|0),-, it may undergo a second-order transition |0); -^ |e); — > |0); 
so as to eventually pick up a phase shift once / is scattered off. 
To make rigorous such considerations, we next prove that the 
reflection coefficients are again, with due replacements, given 
by (O as for setup A. 

ft conserves the total number of excitations. Thus, follow- 
ing a standard approach i27i,i28il . we seek one-excitation sta- 
tionary states of the form 

l*I'ar> = ^ I (ixi/rQ.j(x)c^',(x)|vac>|aiQ'2)i2 

-l-ei|vac)|eQ'2)i2+£2|vac)|Q'ie)i2, (11) 

where 4'a±{x) have a form analogous to Eqs. (|4|i and Q 
thus being specified by parameters {r,ai,a2,b\,b2], (e,) are 
excited-state amplitudes and |vac) is the field vacuum state. 
Thus, for given a '^a) is specified by the 7 complex am- 
plitudes {r,ai,a2,bi,b2,B\,B2] and obeys the SE H'^a) - 
TJk\^a)- Projecting this onto c_[_(x)|vac)|Q'ia'2) gives (we 
henceforth omit subscript a) 

+ iviff'j_{x) + J y 6afiSi5{x-Xi)-vk\p±{x) , (12) 

/=1.2 

where J— V2y is the rate associated with the transition |0) <-^ 
\e). Further projection of the SE onto |vac)|eQ'2) and |vac)|Q'ie) 
immediately yields that for each / = 1, 2 s,- = JSa,Qilf{xi)l(vk- 
cl)q). Replacing these in (fT2l i we are left with {i/'(x), tf/±(x)} 
only 

+ iv4r'j^{x)+J^/{vk-cao)2_j 6afi>]j(xt)6{x-xc)-vk[]/±{x) .{13) 

Subtracting now the equation for i/r- from the i/'+'s one 
yields —iif/'ix) - A:[i/^+(x)-i/'_(x)], which trivially entails that 
-/Ai^'lvf - A:(Ai/'+|rf - Ai/r-lvf) holds as well for each £ - 
1,2. Each Ai^ilvf on the rh.s. of the above can be evalu- 
ated by integrating ( fT3] l over an infinitesimal interval across 
X - xi {( - 1,2), which straightforwardly yields Ai/r±|i^ = 
+ii/v)J^ /(vk-ajo)Sa,Q4fixr)- Thereby 



A^'l 



V vk-ojQ 



6,fi>^{xe) (^=1,2) 



This is identical to Eq. ^ once we set 



vk—ojQ 



(14) 



(15) 



and note that, due to the parabolic dispersion law, in setup A 
m — k/u. 



As (/'(x), besides (fT4] i. must fulfill conditions analogous to 
Eqs. ^ due to the common geometry of setups A and B, we 
conclude that amplitudes {ra} for setting B are identical to dSJ 
with the effective mass and potential height given by k/v and 
(flSl l. respectively. In passing, note that \ra\'^ - 1 showing that 
if/ is absorbed it will be re-emitted with certainty (each atom 
behaves as an effective qubit encoded in {\go}, \gi)})- It is now 
evident that setup B can be used as an emulator of A, thereby 
allowing for occurrence of the CZ gate [cf. Eq. (fTOl iI under 
conditions ( fTSI l. Interestingly, the requirement y = F/k » 1 
now in fact becomes the resonance condition (RC) vk ^ ojo- 
This agrees with luTll . where it was shown that the reflectivity 
of an atomic scatterer becomes unitary in this limit. Also, al- 
though of a different nature the first two requirements in (fT6] l 
are RCs either: The CZ gate thus stems from a combination 
of RCs. As anticipated, setup B can be experimentally imple- 
mented in several different ways, including photonic-crystal 
waveguides with defect cavities II18I1 . semiconducting (dia- 
mond) nanowires with embedded QDs (nitrogen vacancies) 
|19,|20], optical or hollow-core fibers interacting with atoms 
1 2l|] and microwave (MW) transmission lines coupled to su- 
perconducting qubits 1I22I1 . In particular, highly efficient cou- 
pling between the dot and the fundamental waveguide mode 
has been recently achieved in tapered InP nanowires with em- 
bedded InAsP QDs [201. Interestingly, it was recently shown 
II26I1 that a hard-wall BC analogous to ours [cf. Fig. 2(b)] can 
benefit MW photodetection. It is also worth mentioning that 
some features of our scheme are reminiscent of Refs. 12911 and 
Mi- 
Gate working condition. In practice, the incoming / has 
a narrow but finite uncertainty Ak around a carrier wave vec- 
tor ko fulfilling (fTST l. As R is fc-dependent, we assessed its 
resilience against deviations from fco by using the process fi- 
delity F as a figure of merit 11 3111 . In a representative case, 
Ak<5%ko yields F>95%. Also, A^ affects the gate duration 
At according to At~ l/(vAk) ll3lll . which entails a minimum 
time AT,j,in ~ 10 {vk^y^ . If T^ is the system's decoherence time 
the gate works reliably when Atj^w <^ T^, hence the working 
condition reads T^ » 10(vA;o) '. This is matched in realistic 
instances JSlll . 

Conclusions. We have shown a strategy to quasi- 
deterministic ally carry out multi-qubit gates between static 
qubits through single flying buses scattering from them. This 
is effective without demanding post-selection of any kind or 
iteration. The possibility to naturally implement a universal 
CZ gate has been proven in two different setups including ID 
photonic waveguides coupled to atom-like qubits. We believe 
this work can set a significant milestone for future advance- 
ments in the area of distributed QIP as well as in the emerging 
field of quantum optics in ID waveguides. The design of 
a full quantum computing architecture, implementing single- 
qubit operations besides the proposed multi-qubit gates, will 
be the subject of a future comprehensive work 113211 . 
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SUPPLEMENTARY MATERIAL 

The aim of this Supplementary Material is to carry out a de- 
tailed analysis of our scheme's performances by relaxing the 
assumption (made in the main text) to deal with a monochro- 
matic incoming wavepacket. We first quantify the scheme's 
resilience to deviations of the wave vector from the optimal 
value yielding a CZ gate. Specifically, we evaluate how dif- 
ferent is the implemented gate from the ideal CZ by using a 
monochromatic wave packet whose associated wave vector is 
arbitrarily chosen. This can also be regarded as the fidelity 
corresponding to a generic harmonic secondary wave when a 
finite-width wavepacket is sent. Next, under the latter (real- 
istic) hypothesis of dealing with a wavepacket, we work out 
the characteristic time taken by the scattering process. This 
turns out to be a function of the group velocity and wavepacket 
width. Finally, we use jointly the above findings in order to 
define the gate's working condition in terms of the system's 
decoherence time. 



A. Gate fidelity 

Here, we analyze the resilience of the two-qubit CZ gate 
proposed in the main text (setups A and B) to an imperfect 
setting of k, i.e. the wave vector of the flying particle. For set 
X21, X32, « and n' and in the limit y » 1, the optimal wave 
vector ko fulfills [cfr. Eq. (9) in the main text] 



kQX2i = nn, koX32 = (n' + l/2)7T . 



(16) 



The corresponding reflection operator R^^ coincides with the 
ideal CZ gate and can be written as 



^-l-o = |0)l<0|l2 + |l)l<l|Z2 



(17) 



where 1,, X,, Y, and Z, are respectively the identity operator 
and the usual Pauli operators of the /th static qubit. Using 
Eq. (8) in the main text in the limit y » 1, the gate corre- 
sponding to an arbitrary k is given by 

^^ = |0>i<0| I2 + |1>1<1| (e-*"=|0>i0|+e2*'-^|l>2(l|) , (18) 

which in general does not coincide with ( fTTl i. To quantify 
how close is Rk to R^o, we use the process fidelity F ll33ll as 
a figure of merit. This is defined as F - TrOtoA')' where the 
16x16 matrix ;^' fully specifies R^ according to 



^k P^k^Z^ Xmn A„, p A„ 



(19) 



where p is a generic two-qubit state while jA,,,) is a basis for 
operators acting on p [xo specifies Rk^j in full analogy with 
(fT9]l1. We choose each A,„ to be the tensor product between 
two operators taken from the set {t,X,Y,Z]. Hence, through 
the replacements |0X0| = (l+Z)/2 and |1X1| = (ll-Z)/2 in 
Eqs. ( fTTl ) and ( fTSl l one can exactly derive both xo and Xk- 
Upon trace of their matrix product, the gate fidelity can be 
worked out in the compact analytical form 

3 -H 2 cos(2fcx2) - cos(2A;ji:32 ) - 2 cos(2A;a:3 ) 
F^ . (20) 




1.00 



FIG. 3. (Color online) (a) Gate fidelity F against kikf). (b) Zoomed 
view of panel (a) showing that a discrepancy between k and ^o up 
to over 5% can be tolerated in that F stays above the 95% threshold 
(red dashed line in the plots). 



We will next show that the above features of the incoming 
wave packet guarantee that the characteristic time At is the 
same as the one associated with the free-propagating particle, 
i.e. in the absence of the static qubits and mirror. The phys- 
ical mechanism behind this behavior is essentially the one in 
Ref. 113411 . where however a different Hamiltonian model was 
addressed. The main features of the demonstration are basi- 
cally common to both setup A and B. We will carry out the 
demonstration in detail for the latter since this is somewhat 
more involved owing to the second-quantization formalism of 
the associated model. 



In order to assess the gate resilience, it is convenient to 
study F as given in ( l20t as a function of the dimensionless 
parameter A:/A:o, where k^ fulfills ( fTSI l. Thereby, F - \ for 
klk(j = 1 . In Fig. 1 (a), we set n-\ and «' = [cfr. Eq. ( fT6b l and 
plot F against klko- For |A:/^ol S 0.4, F decays as the discrep- 
ancy between k and ^o grows (regardless of the sign). Impor- 
tantly, as can be seen from Fig. 1(b), relative deviations of k 
from the optimal value as large as ~5% still ensure the fidelity 
to stay above 95%, which amounts to a quasi-deterministic 
gate operation. 



B. Gate characteristic time 

Here, we derive the characteristic time At taken by the CZ 
gate proposed in the main text (setups A and B) in order to 
be implemented through a single scattering process. To this 
aim, we will first relax the assumption to deal with a perfectly 
monochromatic incoming wavepacket for the flying particle 
(as such, this would clearly entail an infinite At). We in- 
stead consider a naiTow Gaussian wave packet impinging on 
the static qubits and work out the associated energy-time un- 
certainty principle AE, which is related to the characteristic 
time by the energy-time uncertainty principle A/JAt ~ 1 (we 
recall that we set S = 1). 

As incoming flying-particle wave function ip{x) we take the 
right-propagating Gaussian function 



ip{x)-- 



1 



^/2^Ak 



(21) 



specified by the center xq, the carrier wave vector A:o > and 
the wave-vector width AA:. We assume that xq and Ak fulfill 



|-*:i-xo| > 3Ax , k^)>3Ak , 



(22) 



where Ax is the wave-packet width in position space. The 
first of the above conditions specifies that the incoming wave 
packet is initially fully out of the scattering region (we re- 
call that this lies at xi < x < x^, see Fig. 2 in the main 
text). The second condition sets an upper bound on the wave- 
vector width (~0.3 ko) so as to ensure that the incoming wave 
packet (p{x) comprises only positive-A: secondary plane waves 
(i.e. right-propagating ones). Owing to the latter feature the 
i^(x)'s Fourier transform ip(k) - l/y/lK j dxe^''''^ip(x) thus 
obeys 



Initial state 



If the static qubits are in the state lor) 12 (see main text) and 
the incoming flying particle is described by ( |2TI ). the initial 
system's state reads 



\% 



•J —o 



dx(f(x)c'(x)\wsic)\a)i2 - 



(24) 



We first show that due to (l23l l in jTM c'(x) can be replaced 
with c+(x) (bosonic operator creating a right-propagating pho- 
ton at X, see main text). As c'(x) = cl{x) + c'_{x) we need to 
prove that J ^ dx tf(x)c'_(x)\vsic) ^ 0. By replacing in the latter 

expression the Fourier decomposition c_(x) = J ^ dk e c^ (q 
is a standard photonic annihilation operator in the A:-space) we 
end up with 



f 

•J — Cx 



dx if{x)c'_{x)\ysic) ■ 



<J —oc 

-f 

•J —cc 



dk 



jZ,dxe "^'(fiix) 



|vac) 



dkip{k)c[\v?iC)^0 , (25) 



where in the last step we made use of ( |23] l. Thereby, the initial 
state is well approximated as 



1%) 



kJ — cx 



dx (p{x)c\ {x) I vac) I a) \ 



(26) 



^(fc < 0) ^ . 



(23) 



D. Spectral decomposition 

To eventually calculate AE we need to derive the spectral 
decomposition of I^Pq) with respect to the system's stationary 
states \^ka) (it is now convenient to explicitly indicate the 
dependance of these on k). Use of Eq. (11) in the main text 
yields 

<*I'^Qrl*I'o> = ^ I rfxc/x>^Q,^(x)(^(x)<vac|cj(x)cl(x')|vac) . 
d=± ^-°° 

(27) 

Next, using that the vacuum state vanishes whenever 
any annihilation operator is applied to it we obtain 
(vac|Q(x)c|(x')|vac) = [q(x), c|(x')]. Such commutator is 
zero for d - - because any right-propagating bosonic op- 
erator commutes with any left-propagating one. We thus 



need to calculate only [c+(x), c+(x')]. To this aim, we use 
c+{x) -J dk e''^^Ck and [c/., q/] = 6{k—k') and obtain 



[c^{x),c\{x')]^5{x-x 



27r J_„ 



dke 



ik(x-x') 



6{x—x') i 



(28) 



x' —X 



2 2ti{x-x') 

Hence, the scalar product in (l27l i becomes 

(29) 
The principal value of the improper integral on the r h.s. can be 
worked out in a standard way with the help of the Cauchy in- 
tegral theorem and formula, the Jordan lemma and (l23T l. Note 
that owing to the latter only positive-A: plane waves appear in 
the Fourier expansion of (f{x) whose extension to the com- 
plex plane (f{x + iy) thus vanishes for y — > oo. Hence, the inte- 
gral is well approximated as J^^ dx'(f{x')/(x'—x)^imp{x) and 
thereby (|29] | is simply given by 

dxil,la^(x)ip(x) . (30) 

CO 

To evaluate the r.h.s. of the above equation, we use Eq. (4) in 
the main text along with the first of (l22l i. Owing to the latter, 
ip{x) is negligible for x > xi. It should now be clear that by 
setting jci = the contributions to the integral proportional to 
6{x)-6{x-X2) and 6{x-X2) [cfr. Eq. (4)] can be neglected while 
the first Heaviside step function 6{-x) can be omitted for all 
practical purposes. Thereby, we are finally left with 



mal^o> = 



i_ r 

7.7T J-K 



V2^ J- 



dxe-'"ip(x) = ip{k) 



(31) 



E. Characteristic time 



By recalling that the energy of l^ka) is E = vk the usual 
decomposition in the steady-state basis of the evolved state at 
a time t reads 



mm- 



■■!■ 



dk{^,a\'i'o}e-'''"\^ka)- 



'-fdkif{k)e-''''"\^'ka) 



(32) 

where we have used the result i3T[ . Eq. ( |32] | is analogous to 
the case where the flying particle propagates freely but for the 
replacement of each plane wave with the steady state of the 
full scattering problem corresponding to the same k. Hence, 
while the dynamics will be obviously different the energy un- 
certainty is necessarily given by AE = uAk as in the free- 



particle case. Using AEAt ~ 1, it follows that the charac- 
teristic time is also the same in the two cases and simply reads 



At- 



1 
uAk 



(33) 



These findings are in Une with the content of Ref. 134 

As for setup A (see related section in the main text), the 
demonstration proceeds analogously but is more straightfor- 
ward because of the first-quantization picture. Indeed, one 
assumes an incident wave packet of the form ( |2TI ) fulfilling 
conditions (l22T i. Then, using these, Eq. (|30] | is immediately 
worked out (without the need for dealing with commutation 
rules) and next Eq. dJTT i is found. Eq. ( |32] | still holds aside 
from the replacement of the linear dispersion relation with 
E - k^/{2m) (see main text). The corresponding free-particle 
energy uncertainty now becomes AE - VkAk, where v - k/m 
is the group velocity corresponding to k. Finally, we end up 
with Ar~l/(L'i.AA). 



F. Gate working condition 

To work out the gate optimal working condition, one needs 
to take into account Eqs. ( [2Q] i and ( [33T l jointly. As shown in 
the first section of this Supplementary Material, in a represen- 
tative case the gate fidelity F stays above 95% for Ak up to 
~5%ko. Using ( |33] ). this immediately entails that the shortest 
time taken by the scattering process to reliably implement the 
CZ gate is given by 



ATrr 



1 



v5%kQ 



■lOivko)- 



(34) 



Clearly, in order for the gate to be effective one needs the pro- 
cess to be short enough compared to the decoherence time 
Td (the latter is the typical time to elapse before the system's 
coherence is significantly spoiled by environmental interac- 
tions). We thereby demand that At «: T^. Since this must 
occur under the constraint At > Armin we conclude that the 
gate can work provided that Armin "^ T^, which in the ligtht of 
OH) gives 



r,, » lO(yfco)" 



(35) 



As representative numerical instances, let us consider a GaAs 
nanowire with embedded InAs quantum dots 113511 . By esti- 
mating the group velocity as f = c/ncaAs, where «GaAs - 3.4 
is the refractive index of GaAs and the typical carrier wave- 
length as Ao = In/ko ~ 900 nm the above condition yields 
T^i s> 1.6 ■ 10 '"^s. This is fulfilled since typical decoher- 
ence times of quantum dots of this sort are of the order of 
ps. Next, we address a diamond nanowire where Nytrogen- 
vacancy (NV) centers work as atom-like systems IBol . In such 
a setting, the refractive index is ^2.4 while Aq ^ 640 nm giv- 
ing rj»8-10"'^s. Again, this is very well matched in the light 
of decoherence times of NV centers, which typically exceed 
tens of ns IBTIl. 
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